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Abstract 

We investigate the chiral phase transition at finite temperature (T) in colour SU(iV c = 3) Quantum Chromodynamics (QCD) with 
a variable number of fermions Nf in the fundamental representation by using lattice QCD. For Nf — 6 we study the approach to 
asymptotic scaling by considering lattices with several temporal extensions N,. We then extract the dimensionless ratio T c /Al (Al = 
Lattice Lambda-parameter) for Nf = 6 and Nf = 8, the latter relying on our earlier results. Further, we collect the (pseudo) critical 
couplings /3 L C for the chiral phase transition at Nf — (quenched), and Nf = 4 at a fixed N, = 6. The results are consistent with 
enhanced fermionic screening at larger Nf. The ratio T c /Al depends very mildly on Nf in the Nf — - 4 region, starts increasing 
at Nf = 6, and becomes significantly larger at Nf = 8, close to the edge of the conformal window. We discuss interpretations of 
these results as well as their possible interrelation with preconformal dynamics in the light of a functional renormalization group 
analysis. 
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^ (— | INTRODUCTION Conformal invariance is anticipated to 
l — 1 emerge in asymptotically free non-Abelian gauge theories when 
[the number of flavours exceeds a critical value Nf = Nf . The 
^ approach to conformality for Nf < Nf, close to the edge of the 
conformal window, is in principle associated with a walking be- 
l/~) haviour of the running coupling, which has been advocated as a 
7—1 [basis for strongly interacting mechanisms of electroweak sym- 
^**\ -merry breaking . 

Recent lattice studies ^ focused on the computation of Nf 
,— ' and the analysis of the conformal window itself, either with fun- 

, damental fermions yl or other representa- 

• • tions 

GHEE 

11411 . Among the many interesting results with 
. ^ [fundamental fermions, we single out the observation that QCD 
-with three colours and eight flavours is still in the hadronic 
?— i [phase @], while Nf = 12 seems to be close to the critical 
. 5^ -number of flavours, with some groups favouring conformality 
fllllElzL an d others chiral symmetry breaking U. 

In comparison, much less effort has been devoted to the 
analysis of the phenomenologically relevant subcritical region 



EHQ3Q1]- Here, we would like to learn where and how 



QCD at the edge of the conformal window displays walking, 
and its associated manifestations like separation of scales and 
approximate scale invariance. This Letter is one step in this 
direction. 

Building on the above mentioned results, we decided to con- 
centrate ourselves on Nf < 8, so to be safely in the hadronic 
region, but not too far from the edge of the conformal window. 
A recent study Hfl noted an enhancement of the zero tempera- 
ture ratio (ifnJf)/F 3 , where F is the pseudoscalar decay constant. 
This suggests that Nf = 6 might indeed be the onset of new 
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strong dynamics. 

In this letter we study the thermal transition of QCD with Nf 
= 0,4,6, and combine our findings with those of our early work 
for Nf = 8 JH. We confirm the expected enhanced screening 
when the number of flavours increases, and discuss the interre- 
lation of our results with a possible emergence of a new, pre- 
conformal dynamics. 

As a general remark, we note that using the thermal transition 
as a tool for investigating preconformal dynamics was largely 
inspired by a renormalization group analysis 11911 . as we will 
review below. Further reasons of interest for finite temperature 
studies of large Nf include a connection between the quark- 
gluon plasma phase and the cold conformal region, which might 
lend support to analyses of quark-gluon plasma based on the 
AdS/CFT correspondence. 

BRIEF OVERVIEW OF ANALYTIC RESULTS A second 
zero of the two-loop beta-function of a non-Abelian gauge the- 
ory implies, at least perturbatively, the appearance of an in- 
frared fixed point (IRFP) and the restoration of conformal sym- 
metry |2(J 21]. In colour SU(3) gauge theory with Nf massless 
fermions in the fundamental representation, the second zero ap- 
pears for Nf > 8.05, before the loss of asymptotic freedom 
(LAF)atA^ AF = 16.5. 

One expects that conformality should emerge when the 
renormalized coupling at the would be IRFP is not strong 
enough to break chiral symmetry. This condition provides the 
lower bound Nj- of a so called conformal window in Nf. Ana- 
lytic studies based on the Schwinger-Dyson equation with rain- 
bow resummations |22,|23j|24! or the functional renormaliza- 
tion group method 1 19(1 suggest ~ 12. An all-order pertur- 
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bative beta-function 12511 inspired by the NSVZ beta-function 



of SQCD 12611 has been conjectured, leading to a bound > 
8.25; N c f has also been estimated for different fermion repre- 
sentations 12711 . In addition to the lower bound of the confor- 
mal window, walking dynamics in the preconformal region is 
another interesting subject in relation to proposed scenarios of 
walking technicolor. Instanton studies at large Nf l28Tl claimed 
a qualitative change of behaviour at Nf = 6. 

All the above phenomena happen well into the strong cou- 
pling regime, rendering a perturbative prediction unreliable and 
a non-perturbative analysis mandatory. The genuinely non- 
perturbative lattice formulation of gauge theories is thus a nat- 
ural candidate for this study. 

Recently, the functional renormalization group (FRG) 
method has been applied to finite T QCD with varying num- 
ber of flavours, and the critical temperature for the chiral phase 
transition was obtained as a function of Nf H]. In this T- A?) 
phase diagram, the onset of the conformal window has been 
estimated by locating the vanishing critical temperature. 

Most interestingly, the critical exponents associated with the 
behaviour of the beta function at the IRFP manifest themselves 
also in the shape of the thermal critical line in the vicinity of 
the critical number of flavours N c f . In more detail, the line is 
almost linear with Nf for small Nf, and displays a singular be- 
haviour when approaching Nf. As emphasised by the authors 
m Ref. JH, the result clearly elucidates the universality of the 
critical behaviour at zero and non-zero temperature in the vicin- 
ity of N^. It thus seems a promising direction to extend the 
knowledge of finite T lattice QCD to the larger Nf region, by 
using the FRG results as analytic guidance. 

In this work we investigate the thermal chiral phase transi- 
tion for Nf = 0, 4, 6, 8 colour SU(N C = 3) QCD by using lat- 
tice QCD Monte Carlo simulations with staggered fermions. 
Nf = 6 is expected to be in the important regime as suggested 
by the results in Refs. 151 12^1 . This work includes the first study 
of Nf = 6 staggered fermions at finite T, and it provides an 
important ingredient to a broader project that studies the emer- 
gence of the conformal window in the T -Nf phase diagram. In 
addition to Nf = 6, we compute the (pseudo) critical coupling 
for Nf = (quenched) and Nf — 4 at N, = 6, and use the results 
from Ref. J5[ for Nf = 8. 

In short, this work explores a largely uncharted territory: the 
chiral transition of strong interactions at high temperature, and 
large number of light flavours. Our goal is to observe, and 
understand possible qualitative differences with the very well 
known behaviour of QCD thermodynamics. To this end, we 
have for the first time collected and analysed results for differ- 
ent number of flavours, and the same lattice action, so to be able 
to meaningfully compare them. We ask the question: are we 
still finding just small differences among theories with different 
number of flavours, or are we going to observe some significant 
trend? As we will see, our exploratory, and in many respects 
qualitative analysis, will indicate that Nf = 6, and even more 
Nf = 8, are serious candidates for a different chiral dynamics. 
Obviously, our observations call for detailed quantitative stud- 
ies which are already underway. 

SETUP Simulations have been performed by utilising the 



publicly available MILC code 12911 . The setup explained below 
is the same as the one used for Nf = 8 in Ref. |5J. We use 
an improved version of the staggere d action, the Asq tad action, 
with a one-loop Symanzik 13Cj bill and tadpole 13411 improved 
gauge action, 

S = -^TrlogM[am, U,u ] + J] MslW ~ Uq] , (1) 

i=P,r,pg 

where g^ is the lattice bare coupling, and /?,■ are defined as 
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The plaquette coupling /3 P = lO/g^ = Pl is a simulation in- 
put. The M[am, U, Uq] in Eq. (Q]) denotes the matrix for a single 
flavour Asqtad fermion with bare lattice mass am, and Uq rep- 
resents the trace of the ordered product of link variables along 
C,, for the lxl plaquettes (i = p), the 1x2 and 2x1 rectangles 
(i = r), and the lxlxl parallelograms (i = pg), respectively 
- all divided by the number of colours. The tadpole factor uq is 
determined by performing zero temperature simulations on the 
12 4 lattice, and used as an input for finite temperature simula- 
tions. 

To generate configurations with mass degenerate dynamical 
flavours, we have used the rational hybrid Monte Carlo algo- 
rithm (RHMC) [32], which allows to simulate an arbitrary num- 
ber of flavours through varying the number of pseudo-fermions. 
Simulations for Nf — 6 have been performed by using two 
pseudo-fermions, and subsets of trajectories for the chiral con- 
densates and Polyakov loop have been compared with those ob- 
tained by using three pseudo-fermions with the same Monte 
Carlo time step dr and total time length r of a single trajec- 
tory. We have observed very good agreement between the two 
cases for both evolution and thermalization. We have moni- 
tored the Metropolis acceptance and reject ratio, and adjusted 
t = 0.2- 0.24 and dr = 0.008 - 0.018 to realize the best perfor- 
mance. For each parameter set, we have collected a number of 
trajectories ranging from a one thousand to five thousand - the 
latter closer to the critical region. 

The focus of this letter is the analysis of the chiral transi- 
tion. The fundamental observable is then the order parameter 
for chiral symmetry, the chiral condensate: 
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where N s (N t ) represents the number of lattice sites in the spa- 
tial (temporal) direction and is the temporal link variable. 
We have also measured connected and disconnected chiral sus- 
ceptibilities, 
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and we have considered the logarithmic derivative, 
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and 
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As discussed in previous work ||5, 33], R„ is a probe of chiral 
symmetry which is particularly useful for numerical investiga- 
tions. To further characterize the critical region, we also mea- 
sured the Polyakov loop, 



L = 



N c Ni 
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and tr c denotes the trace in colour space. 

The temperature T is related to the inverse of the lattice tem- 
poral extension, 



1 



a(J3 L ) ■ N, 



(10) 



We measure (ijnfr), the chiral susceptibilities and L at various 
temperatures. The output of this measurement is the (pseudo) 
critical coupling y6 L ° for the chiral phase transition for a given 
value of N t . We underscore that all the measurements of the 
pseudo-critical couplings which will be used in our discussion 
are only based on fermionic observables. 

RESULTS All results have been obtained for a fermion bare 
lattice mass am = 0.02. In Figs. Q] and [2] the expectation val- 
ues of the chiral condensate a 3 (t//i/f), and the Polyakov loop 
L are displayed as a function of /3l for several N,, respec- 
tively.lt is found that different give a different behaviour 
of a 3 (ifnf/) and L. In particular, this indicates that their rapid 
crossover with increasing /?l is not to be attributed to a bulk 
transition. The asymptotic scaling analysis below will confirm 
that it corresponds instead to a thermal chiral phase transition 
(or crossover) in the continuum limit. Q 

For N t = 4, it is possible to extract /? L C = 4.65(25) from 
the peak position of the first derivative - a'd^ili) Idfi^. We 
note that at this level of accuracy we cannot disentangle the 



A note on the mass dependence is in order. According to the Pisarski- 
Wilczek scenario, the most likely possibility for Nt >. 3 is a first order chiral 
transition in the chiral limit. Standard arguments indicate that first order phase 
transitions are robust against explicit breaking: when introducing a bare quark 
mass, then, we expect a first order phase transition which will eventually end in 
a genuine singularity at some critical point (T c , m c ). By further increasing the 
bare mass, the transition will turn into a crossover. In the unexpected situation 
of a second order transition in the chiral limit, any nonzero quark mass will 
immediately produce a crossover. Even in the case of a first order transition, 
though, a finite lattice will turn it into a crossover. All in all, as in any lattice 
study, we are never dealing directly with a genuine criticality. Rather, we are 
locating a pseudo-critical point, and only by considering several masses and 
several volumes, we can assess with confidence the nature of the phase transi- 
tion in the infinite volume and in the chiral limit. Discriminating among these 
different behaviours is however beyond the scope of this study. 



peak position for the derivatives of the chiral condensate and 
the Polyakov loop dL/dp^. 

For N, = 6, we find a small jump between /?l = 5.0 and 5.05, 
where the Polyakov loop also shows a significant enhancement. 

For N, = 8, the Polyakov loop L shows a clear signal as in- 
dicated in Fig. [3] In particular, we observe a drastic increase of 
the Polyakov loop L around 5.2 < /?l < 5.3. The histogram of 
the chiral condensate around the rapid increase of the Polyakov 
loop is shown in Fig. [4] At /?l = 5.2, the histogram exhibits a 
broadening that suggests the increase of fluctuations around the 
pseudo-critical point. 

For N t = 12 the results are particularly smooth. Note that 
in this case the aspect ratio is only two, and larger volumes 
would be required to reach a comparable clarity in the signal, 
although the onset for the Polyakov loop at /? = 5.50(5) is still 
appreciable. 

The results for R„ are collected in Fig|5] R„ has a clear 
jump for/?,. = (4.65(5), 5.05(5), 5.20(5)) for N, = (4,6,8). For 
Af, = 12, the behaviour is again rather smooth. We have then 
searched for an inflection point in R„ by fitting to a hyperbolic 
tangent in several intervals, varying the extrema of the integra- 
tions. The results are reasonably stable in the errors. We will 
quote as a central value the average of the fit results within dif- 
ferent intervals, and for the error, we will use the conservative 
estimate of half the difference between the largest and small- 
est result, obtaining y8 c = 5.46(14). All values of the (pseudo) 
critical lattice coupling y8 L c are summarized in TableQ] 

These results can be analyzed and interpreted in terms of the 
two-loop asymptotic scaling law. Let us consider the two-loop 
lattice beta function, 
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for fundamental fermions in SU(N C ) (C 2 [G], C 2 [F], T[F]) = 
(N c , (N 2 - \)I(2N C ), 1/2). From Eq. <□} we obtain the well 
known two-loop asymptotic scaling law, 
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Here Al is the so called lattice Lambda-parameter, and /3l = 
2N c /g 2 . 

The above relation is valid in the massless limit. In the fol- 
lowing, we will use it to analyze results obtained at finite mass. 
This assumes that the shift of the (pseudo) critical coupling in- 
duced by a non-zero mass is smaller than other errors. This 
assumption should ultimately be tested by performing simula- 
tions with different masses and extrapolating to the chiral limit. 

The equation dTOb can be written as 



1 = aHa l <«). 



(15) 
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The left-hand side is a given number, and we have obtained the 
corresponding /? L C by lattice simulations. Hence, the quantity 
T c /Al on the right-hand side of Eq. ( TOl ) can be extracted, and 
must be unique as long as the asymptotic scaling law Eq. dT4b 
is verified for a given j8 L c . 

If we use the lattice bare coupling to carry out the above pro- 
gram, we notice appreciable scaling violations. Indeed, Eq. (TBI) 
holds true up to non-universal scaling-violating terms. One pos- 
sibility is then to follow earlier work l35ll and to parametrize 
scaling violations as 

A L a(J3 h ) ee [A Lfl (/? L )] 21oops (l + /i[A L flOS L )]L ops ) ■ (16) 

where [A L a(/? L )] 21oops is defined by Eq. (14]). 

Alternatively, we can trade the bare lattice coupling pl for 
the boosted coupling introduced in our previous work |5j] g — 
y/2N c /lO ■ gL- We will show below that this prescription leads 
to a rather accurate two-loop scaling, equivalent to the en- 
hancement of the scaling behaviour obtained by considering 
Eq. ffiM . Our prescription is similar in spirit to the Parisi- 
Lepage-Mackenzie 13411 boosted coupling. Our boosted cou- 
pling, however, cannot be derived in (tadpole improved) per- 
turbation theory. In our future work, we plan to perform zero 
temperature measurements which will allow an independent es- 
timate of the lattice spacing, and a complete discussions of scal- 
ing, and asymptotic scaling. Hopefully, this will shed light 
on the reasons why our simple approach works so well. At 
this stage, admittedly, it remains a heuristic, ad-hoc prescrip- 
tion which effectively incorporates the scaling violations in the 
range we have explored. This warning issued, we will continue 
the discussion by using our boosted coupling. 

In Fig. [6] we show the A^ 1 - A^a(J3^) plot. The slope of 
the line connecting the origin and the data points corresponds 
to T c /Al, The N, = 6, 8, and 12 points have a common slope 
to a very good approximation, while the N t = 4 result falls on a 
smaller slope. 

The latter is interpreted as a scaling violation effect due to 
the use of a too small A^. The existence of a common T c /A^ 
for N, > 6 indicates that the data are consistent with the two- 
loop asymptotic scaling Eq. (fT4l . confirms the thermal nature of 
the transition and that Nj — 6 is outside the conformal window, 
as expected from a previous Nf — 8 study JH. A linear fit 
provides T c /A^ = 1.02(12) x 10 3 , which can be interpreted as 
the value in the continuum limit for Nf — 6 QCD. Note that 
this compares very well with the result obtained by using f3^ 
obtained by N, — 8 simulations. In the following, we can then 
use it as a representative result for six flavours. 

In order to have a more complete overview, we have per- 
formed simulations for the theory with Nf — (quenched) and 
Nf = 4, only at N t = 6. These theories are of course very well 
investigated, however we have not found in the literature results 
for the same action as ours. Table Q] shows a summary of our 
results for the (pseudo) critical coupling /3 L C of the chiral phase 
transition at finite temperature for Nf = 0, 4, 6, and 8 - the 
latter from Ref. H]. 

DISCUSSION In Fig. El we display the (pseudo) critical 
values of the lattice coupling g c = J2NJ]}£ from Table Q] in 
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Figure 1: The chiral condensate a (ifrt/r) for Nf = 6 and am = 0.02 in lattice 
units, as a function of for N t = 4, 6, 8, and 12. Error-bars are smaller than 
symbols. 
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Figure 2: The Polyakov loop L for Nf = 6 and am = 0.02 in lattice units, as a 
function of for N, = 4, 6, 8, and 12. Error-bars are smaller than symbols. 



the Miransky-Yamawaki phase diagram. 

Consider the N, = 6 results: it is expected that an increasing 
number of flavours favors chiral symmetry restoration. Indeed, 
we find that, on a fixed lattice, the (pseudo) critical coupling 
increases with Nf in agreement with early studies and naive 
reasoning. The precise dependence of the (pseudo) critical cou- 
pling on Nf at fixed Af, is not known. It is, however, amusing 
to note that the results seem to be smoothly connected by an al- 
most straight line: the brown line in the plot is a linear fit to the 
data. Comparing the trend for Nf = 6 to the one for Nf = 8, for 
varying N, one can infer a decreasing in magnitude (and small) 
step scaling function, hence a walking behaviour. Further study 
is needed at larger Nf, and by using the same action used for 
Nf = - 8, to confirm or disprove it. 

Next, we study the Nf dependence of the ratio T c /Al and 
related quantities. We recall that the simulations for Nf - 4 and 
Nf = have been performed by using only N t = 6. Hence, in 
these two cases, the results will hold true barring strong scaling 
violations at N, = 6. We note that in a previous lattice study 
with improved staggered fermions I36J, asymptotic scaling was 



Table 1: Summary of the (pseudo) critical lattice couplings f}£ for the theories with Nf = 0, 4, 6, 8, am = 0.02 and varying N t = 4, 6, 8, 12. All results are 
obtained using the same lattice action. 
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Figure 3: Zoom-in of the chiral condensate cr(if/t//) and the Polyakov loop L 
shown in Figs.^and[2]in the critical region at N, = 8, with spatial volume 24 3 . 



Figure 4: Distribution of the chiral condensate a (ipif/) for Nf = 6, am = 0.02 
and spatial volume 24 3 , in the vicinity of the chiral phase transition at N t = 8. 



indeed observed using a boosted coupling for N, > 6 for < 
Nf < 4. 

Ideally, we would like to convert our results to T c /A^ s . Un- 
fortunately, to our knowledge, the conversion from Al to A^g 
for a generic number of flavours is only available for Wilson 
fermions 13711 . 

Here we consider a simplified procedure, aiming at capturing 
at least the basic features induced by setting a UV scale. For this 
purpose, we introduce a reference coupling y8 L let and an associ- 
ated reference energy scale A, e f. Then Eq. ( fT4l is generalized 

as 



the limit 
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At leading order of perturbation theory b[ 
related via 
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-> 0, Al and A re f are 
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This equation would be analogous of the ratio Al/Ams derived 
in l37ll for Wilson fermions up to a further linear dependence 
on Nf in the numerator of the exponent. In a nutshell, the differ- 
ence originates from the fact that we are fixing a bare reference 
coupling j6 L lef , which will be specified later. Notice that by con- 
struction A re f reproduces the lattice Lambda-parameter Al in 



Aref(y8L ef ^0) = A L (l+<9(l/y6 L C )) 



(19) 



In summary, when trading Al for A re f , we are moving towards 
a more UV scale. 

Let us consider first T c /Al. The values of T c /Al are sum- 
marized in Table [2] and plotted in Fig. [8] The ratio does not 
show a significant Nf dependence in the region < Nf < 4, it 
starts increasing at Nf = 6, and undergoes a rapid rise around 
Nf = 8. The chiral phase transition would happen when T be- 
comes comparable to a typical energy scale M x = CAl- The 
nearly constant nature of T c /Al in the region Nf < 4 indicates 
that the role of such energ y sc ale is not significantly changed 
by the variation of Nf (see 13811 for a detailed discussion of this 
point. ) In turn, the increase of T c / Al in the region Nf > 6 might 
well imply that the chiral dynamics becomes different from the 
one for Nf < 4. Indeed, a recent lattice study |[l5Tl indicates that 
Nf — 6 is close to the threshold for preconformal dynamics. 

We now consider T c /A le f. The Nf dependence of the ratio 
R(N f ) = (T c /A Kf )(Nf) is shown for several £ L ref in Fig.|9] where 
the vertical axis is normalized by R(0) = (T c / A K f)(Nf = 0) for 
each y6 L ret . T c /A K f is now a decreasing function of Nf for a 
larger jS L ref , i.e. for a more UV reference scale A re f. This result 
is consistent with the FRG study Il9ll . where the decreasing 
T c (Nf) has been obtained by using the r-lepton mass m T as a 
common UV reference scale with a common coupling a s (m T ). 

The A re f scale associated with a /? L ref » where /3* is eval- 
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Figure 5: The ratio of scalar and pseudo-scalar contributions to the susceptibil- 
ity, defined in Eq. j6) as a function of /?l- 
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Figure 6: The thermal scaling behaviour of the (pseudo) critical lattice coupling 
y3 L c . Data points for Al a(fl£) at a given l/N, are obtained by using /3/ from 
Table [T]as input for extracting Al in the two-loop expression Eq. H4\ . 

The dashed line is a linear fit with zero intercept to the data with N t > 4. 



uated at the infra-red fixed point should provide a UV scale 
well-separated from the IR dynamics. If we assume the lower 
bound of the conformal window to be Nj ^12, the two-loop 
beta-function leads to /?„ = -2N c bi/b - 0.63. Indeed Fig. 
shows that the decreasing nature of (T c / A re f)(A//) is still weak at 
y8 L ref = 1.0. In the limit y8 L ref -> 0, r c /A ref reproduces Fig. [8] and 
the resultant increasing feature should be attributed to the van- 
ishing of Al due to infra-red dynamics. We also notice that y6 L ref 
must always be smaller than j3 at the UV cutoff, /3\jy = /? L C (A//). 
As shown in Table [T] the lowest value of the (pseudo) critical 
coupling is given by p£(N f = 8,N, = 6) = 4.1125 + 0.0125, 
hence we constrain our analyses to y 6 L J ef < 4.0. In summary, 
Figs.|8]and|9]together show the effects of shifting the reference 
scales from the IR to the UV. 

With the use of a UV reference scale, we should observe the 
predicted critical behavior IU9I1 



T c (N f )=K\Nf-N c f \- 1/e . (20) 
By choosing the critical exponent 9 in the range predicted by 



Figure 7: (Pseudo) critical values of the lattice coupling g c = ^j2N c //3^ for 
theories with Nf = 0, 4, 6, 8 and for several values of N, in the Miransky- 
Yamawaki phase diagram. The dashed (brown) line is a linear fit to the N t = 6 
results. 



FRG: 1.1 

N' 



< 1 /\9\ < 2.5, our data are consistent with the values 
- 9(1) for /3 L ref = 4.0 and N c f = 11(2) for j8 L ref = 2. We plan 
to extend and refine this analysis in the future, and here we only 
notice a reasonable qualitative behaviour. 



Table 2: T c /Al for several Nf. Results are obtained by using the same lattice 
action. For Nf = 6, we have used the N, = 8 result as a representative value. 
The values for Nf = 8 are extracted from Ref. @|. 



N f 


TJA L 





600 ± 34 


4 


620 + 28 


6 


1000 + 92 


8 


2098 ± 191 



SUMMARY We have investigated the chiral phase transi- 
tion and its asymptotic scaling for Nf = 6 colour SU(3) QCD 
by using lattice QCD Monte Carlo simulations with improved 
staggered fermions. This study provides an important ingredi- 
ent to a broader project that studies the emergence of a confor- 
mal window in the T - Nf phase diagram. We have determined 
the (pseudo) critical lattice coupling yS L c for several lattice tem- 
poral extensions N,. We have extracted the dimensionless ra- 
tio T c /Al (Al ^Lattice Lambda-parameter) for the theory with 
Nf = 6 using two-loop asymptotic scaling. The analogous re- 
sult for Nf = 8 has been extracted from Ref. J5|. Tc/A^ for 
Nf = and Nf = 4 has been measured at fixed N, = 6, barring 
asymptotic scaling violations. Then we have discussed the Nf 
dependence of the ratios T c /Al and T c /A le f, where A le f is a UV 
reference energy scale, related to Al as in Eq. ( TT8l . 

We have observed that T c /Al shows an increase in the region 
Nf = 6-8, while it is approximately constant in the region Nf < 
4. We have discussed this qualitative change for Nf > 6 and a 
possible relation with a preconformal phase. We repeat that all 
results have been obtained by working at one value of the quark 
mass and this is a potential weakness of our calculations. 



6 



5 

















1 






























4- 


4- 



















2500 
2250 
2000 
1750 
1500 
1250 
1000 
750 
500 
250 




Figure 8: The ratio 7V/Al, for Nf 
am = 0.02. 
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Figure 9: The Nf dependence oiR(N f )/R(0) for several finite fixed /? L ref . Here, 



R(N f ) = (r c /A ref )(W/). The limit ^ 



reproduces the results shown in 



Fig.[8]up to a renormalization factor and and up to a corrections <5(1//J. C ). 



The ratio T c /A m f is a decreasing function of Nf. This be- 
haviour is consistent with the result obtained in the functional 



renormalization group analysis 11911 . where a common UV ref- 
erence scale was used to study the chiral phase boundary in the 
T - Nf phase diagram. 

Next steps of the current project involve a scale setting at zero 
temperature by measuring a common UV observable. It would 
also be desirable to have the relation between Al and A^ s for 
our action. 

This, together with a more extended set of flavour numbers, 
will allow a quantitative analysis of the critical behaviour. We 
expect the resultant T c - Nf phase diagram to play an essential 
role in the study of the conformal window. 
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